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: ABSTRACT 

(N 

O '■ COBE has provided us with a whole-sky map of the CBR anisotropics. However, even 
O _ 

^ . if the noise level is negligible when the four year COBE data are available, the cosmic vari- 
ed 

G^ ■ ance will prevent us from obtaining information about the Gaussian nature of the primordial 



Ph| fluctuations. This important issue is addressed here by studying the angular bispectrum of 

the cosmic microwave background anisotropics. A general form of the angular bispectrum 
is given and the cosmic variance of the angular bispectrum for Gaussian fluctuations is cal- 



^ ■ culated. The advantage of using the angular bispectrum is that one can choose to use the 



multipole moments which minimize the cosmic variance term. The non-Gaussian signals 
in most physically motivated non-Gaussian models are small compared with cosmic vari- 
ance. Unless the amplitudes are large, the non-Gaussian signals are only detectable in the 
COBE data in those models where the angular bispectrum is fiat or increases with increasing 
multipole moment. 

Subject headings: cosmic microwave background — cosmology: theory 



The cosmic microwave background radiation (CBR) provides one of the most useful tools 
for studying the primordial fluctuations that seed large scale structures today. Valuable 
information about the physical processes that generate primordial fluctuations in the early 
universe can be obtained by studying the statistics of temperature anisotropies in the CBR: 
are they Gaussian or non-Gaussian? Cosmic inflation (Guth 1981; Linde 1982; Albrecht 
& Steinhardt 1982; Liddle & Lyth 1993) predicts a Gaussian pattern of anisotropies on 
all angular scales. Spontaneous symmetry breaking, on the other hand, will lead to the 
formation of topological defects, and the pattern of temperature anisotropies produced by 
defect-networks tends to be non-Gaussian (Bouchet et al. 1988; Turok & Spergel 1991; 
Turner et al. 1991). After COBE's detection of CBR anisotropies on large angular scales 
(Smoot et al. 1992), both skewness and the three point temperature correlation function 
(which contain the lowest order deviations from a Gaussian) were proposed to test the 
Gaussian nature of the CBR (Luo & Schramm 1992, 1993; Falk et al. 1993). However, there 
is a problem in applying these statistics to test for Gaussianity in the COBE map even with 
all four years of data available (Hinshaw et al. 1993): namely the cosmic variance (Abbot & 
Wise 1984; Scaramella & Vittorio 1990, 1991; Cayon et al. 1991; White et al. 1993), i.e. the 
fact that one can only measure CBR temperature anisotropies in one single universe, which 
introduces theoretical uncertainty in the skewness and three point function. Theoretical 
studies (Srednicki 1993) and Monte-Carlo simulations (Scaramella & Vittorio 1990, 1991) 
show that the cosmic variance is large. These results cast doubt on whether one could ever 
get any information about the non-Gaussian nature of CBR from large scale experiments 
such as COBE, simply because of the limitations from cosmic variance. The goal of this 
paper is to address this important issue in the context of the angular bispectrum. We will 
introduce and discuss several properties of the angular bispectrum, including the relationship 
to the bispectrum of the perturbations of the gravitational potential, our modeling of the 
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angular bispectrum, and the cosmic variance of the angular bispectrum when the fluctuations 
are Gaussian. Limits on model parameters and the observability of the angular bispectrum 
are also discussed. 



1 Angular Power Spectrum 

Before we move on to discuss the angular bispectrum, let us briefly review two point statistics 
first. In the case when the temperature anisotropics are Gaussian, the two-point statistic 
is all we need to quantify the temperature anisotropy patterns observed in the experiments 
(such as COBE). 

The GBR temperature anisotropy is a 2D random field defined on the two-sphere. One 
can perform a spherical harmonic expansion of the temperature anisotropy on the sky: 

AT, 

^0 Ira 



Y.aTYne.4>). (1) 



Then the spherical harmonic expansion coefficients a™, / ^ 0, are random variables, and 
their statistics (usually called angular statistics) will completely specify the statistics of the 
temperature anisotropy itself. The angular power spectrum Ci, is the Ith component of the 
Legendre polynomial expansion of the two-point temperature correlation function (Bond and 
Efstathiou, 1987): 

Here the angled brackets denote an ensemble average and Ci is related to the spherical 
harmonic expansion coefficients through: 

(a^ap^'*) = aV^n^V (3) 

In realistic experimental settings there are two major factors which modify the theoretical 
angular power spectrum defined in Eq. (^). One is the finite beam effect, which filters all 



the high multiple moments. For large scale experiments such as COBE, which have a finite 
beam and full sky coverage, the angular power spectrum is modified as (Silk & Wilson 1980; 
Bond & Efstathiou 1984): 

Ci = cm = Ci exp(-Z(Z + 1)(7'), (4) 

where a is the Gaussian beam width. The other important factor is that one can only 
measure CBR fiuctuations in a single universe (our own universe), thus the observed Cf^^ 
is still a stochastic quantity, which has intrinsic fiuctuations. Here C°^^ is the rotationally 
invariant sum 

m=—l 

The cosmic variance of Cf^^ is then given by: 

= iicr?) - ((^r))' = ^^i- (6) 



For a Harrison-Zeldovich (H-Z) primordial spectrum, Q is given by (Peebles 1982) 

Ci 



(7) 



in 5/(/ + l)' 

where Q is the rms ensemble-averaged quadruple. For large I, af oc is a rapid decreasing 
function of I. The cosmic variance dominates the lower multipole moments, and for high l- 
multipoles it is negligible. We expect this is also true for the cosmic variance of the skewness 
and three point function, and we would expect that the theoretical uncertainties due to 
cosmic variance will be reduced by subtracting out lower multipole moments. We will return 
to this point in section 3. 

2 Angular Bispectrum 

If the CBR is Gaussian, the odd moments of the spherical harmonic coefficients af* will van- 
ish. Thus, the non- vanishing odd moments will be a clear signature of deviation from Gaus- 



sianity. The lowest order of such moments, the angular bispectrum, B^{limi,l2m2.,hm^), is 
defined as: 

B^{hrm^,m,^^m^) = (ar^a^^). (8) 

The three point temperature correlation is related to the angular bispectrum through: 

AT " AT AT 
C(/,m,n) = ( — (/) — (m) — (n)) = 
-^o -^o -^o 

i?3(/imi,/2m2,/3m3)l^r(0^r(^)C(^)- (9) 

li,mi 

Because the three point function is rotationally invariant, the angular bispectrum B^{limi, (2^2, 
is non-zero only if kmi, i = 1, 2, 3 satisfy the following conditions [|: 

(1) hjhyh satisfy the triangle rule, i.e. k < \lj — 41, 

(2) h + I2 + h = even, and 

(3) mi + 7712 + ^3 = 0. 

The angular bispectrum can be calculated for any given non-Gaussian model. In this 
paper we consider only the case where the perturbations are adiabatic, so that the temper- 
ature anisotropy is related to the gravitational (p at the last scattering surface through the 
Sachs- Wolfe formula (Sachs & Wolfe 1967), 

= t (10) 

T 3 ^ ' 

The three point temperature correlation function is related to the bispectrum of the gravi- 
tational potential 0, P(j,{ki, k2) = (0/ci0fc2 0fc3)(ki + k2 + ka = 0), through 

^t{1, m,n) = ^^.j P,{k„ k2)e^n^i.^^^^M.h^^^^, + k2 + h)^^^^^, (H) 



^ These conditions can be derived by choosing special beam configurations. Condition (1) is a result of 
the rotational invariance of ^ when two beams coincide. Condition (2) results from the invariance of ^ under 
spatial inversion: I — )■ — ^, rh —rh, h —h. Condition (3) results from choosing beam configurations so 
that the rotational axis is perpendicular to the plane defined by the three beam directions. 
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where t]q = 2c/ Hq is the distance to the last scattering surface and l,7h,n are the beam 
directions. The general relation between the angular bispectrum and the bispectrum 
is rather complicated and we will not give it here. The expression simplifies if P^{ki,k2) 
depends only upon the amplitudes of ki and k2. This is what we expect for the angular 
bispectrum from non-standard infiationary scenario after proper symmetrization. In this 
case, by using the following expansions, 



A7iY.i'ji{kr^,)Y.Yum: 



and 



we find that 



277r2 



dkidk2klklP^{ki, k2)jl{ki)jl{k2)A{lil2l3,mim2 - mg). 



(12) 



(13) 



(14) 



where 

A{lil2h-,rnim2 - 



(2/i + l)(2/2 + l) 



1/2 



Cihhh, 000)C(/i/2/3, mim2 - m^) (15) 



47r(2/3 + 1) 

and C{lil2l3,mim2m3) are Clebsh-Gordan coefficients. The expression 

Alhhh, Tnim2 — m3) is nonzero only if mi+m2+m3 = 0, h + h + h = even and /i, I2, h satisfy 
the triangle rule. Thus, the rotational invariance of the three point function is guaranteed. 

Let us recall that the angular power spectrum Ci is related to the power spectrum P^{k) 
of fiuctuations in through (Kolb & Turner 1990) 

2 



Ci = -^l dkeP^{k)jf{kvo). 



(16) 



The functional form of B^ given by Eq.(|T3p and Eq.(|TB]) motived us to model the angular 
bispectrum as the following: 



Bs{limi, 12^2,131713) = C3/2(0)v4(/i/2/3,mim2 - ma) x 
{a[Cih)Cil2) + Cil2)Cil3) + Cil3)Cih)] + P[Cih)Cil2)Cil3)y}, 



(17) 



where a,/?, 7 are three dimensionless parameters, C{1) = Ci/C{0) is the normahzed angular 
power spectrum and C(0) = ^EG(2/ + 1). For COBE, where the FWHM beam width is 
7°, C(0) = 4.63Q2 and C{1) = 

Apart from guaranteeing the rotational invariance of the three point function, the bispec- 
trum modeled by Eq.(p!7D has another advantage: several physically motivated non-Gaussian 
scenarios give distinctive predictions for the constants a, (3, 7. These constants can be calcu- 
lated directly from Eq.(]I^ and Eq.(|T6|) once the perturbation bispectrum P^{ki, ^2) is given 
for a specific non-Gaussian model. In practice, it was found that in several non-Gaussian 
models of cosmological interest, it is easier to derive an expression for the three point corre- 
lation function ^(fi, f2, fs) in terms of the two point function C{d) than to find an expression 
for P(j,{ki, (Luo 1993). It is convenient to use the normalized two point functions ip and 
three point function 77, where 

nh-h\) = C2{h.h)/C2{Q). ^(o) = i, (18) 

^(ri, f2, h) = i{h, h)/C2{Qf-\ r/(0) = /is, (19) 

and /X3 is the skewness. The theoretical predictions for three point function in specific 
scenarios are the following: 
(1) Inflation 

Various non-standard infiation models will generate a non-zero three point correlation 
function (Falk et al. 1993; Luo 1993; Gangui et al. 1993). The generic form of the three 
point function in most infiationary models is 

V{ri,r2,r3) = y (^(|n-''2|)^(|ri-r3|)+^(|^i-^2|)^(|r2-r3|)+^(|fi--f3|)?/'(|f2-5^^ (20) 

where /i3 is a dimensionless constant. For slow- roll inflation models with one field and a 
cubic self-coupling, /X3 ~ 10^^ ( Falk et al. 1993). Non-linear gravitational evolution will 
produce a three point function of similar form with fi^, ~ 0.01 (Luo & Schramm 1993). 
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(2) xl fields 

Consider the Xn field Y — Yli=iXf, which describes the 0{N)a model of global topo- 
logical defects (Turok & Spergel 1991) in the large N limit, with n — AN. The three point 
function is found to be 

vih, h, h) = ^Jm{'^{\h - r2\m\h - rM\r2 - h\)f^- (21) 

(3) Late-time Phase Transitions 

In this model (Hill et al. 1989), due to the conformal invariance of the system at the crit- 
ical point (Polyakov 1970), the three point function has the following form (Luo & Schramm 
1993): 

^(ri, f2, h) = "(-^dri - f2|)'0(|'^i - ^3|)'0(|^2 - h\)f, (22) 

where a is a dimensionless constant of order unity. 

Once the three point function is found, the angular bispectrum and constants a, /3, 7 
can be found from Eq. (9) by expanding the three point function and two point function in 
spherical harmonics. Non-standard inflation and non-linear gravitational evolution predict 
/3 = 0, q; = /is/S, where /X3 is the skewness; Late-time phase transitions (LTPT) will predict 
q; ~ 0, 7 ~ 1 and /? of order unity; 0{N) a models predict a ~ 0, 7 ~ 1/2 and /? ~ .JqJn 
(Luo 1993). Another possibility is that the underlying CBR fluctuations are non-Gaussian 
only on smaller angular scales. In this case, the angular bispectrum will increase slowly 
with the increasing of I, and will peak around Ic corresponding to the characteristic scale 
9c where the CBR is highly non-Gaussian. If 9c ^ 1°, then Ic ~ 100. Since the observed 
anisotropics in large angle experiments (such as COBE) are the beam-smoothed fluctuations, 
we would expect the angular bispectrum to be fiat, i.e. a = 0, 7 = 0, at low /. Thus, if future 
analysis of COBE map reveals a non- vanishing angular bispectrum, fitting the experimental 



8 



results to the functional form given by Eq.(]T7D will help to discriminate between different 
non-Gaussian models. 



3 Cosmic Variance of the Angular Bispectrum 

The GBR anisotropies are stochastic in nature, which will give rise to a theoretical un- 
certainty in the particular realization of our sky. This effect is severe in large scale GBR 
experiments, and limits our ability to extract information about non-Gaussian nature of 
GBR from current experimental data. In this section we will address this important issue 
and calculate the cosmic variance of the bispectrum when the fluctuations are Gaussian. We 
will show how to minimize the cosmic variance by choosing the appropriate pairs of multipole 
moments. In order for the non-Gaussian signals to rise above the cosmic variance, bounds 
on the model parameters 7 of the bispectrum defined in the previous section are also 
discussed. 

The cosmic variance of the angular bispectrum is given by 

4 = {a-^aT^al^aTraTraTr) (23) 
When ap is Gaussian, by using Wick's theorem, the expression above reduces to 

C3 = Clj^Cl^Cl.^ + 85l-^^l2l^5mim2m3,oCf^ + (^^-^ij ["^mi ,-m2 + '^mim2 ] C'ii C'jg + 

^I2h[^rn2,-m3 + Sm2m3]C'^2^h ~^ ^hhi^rrii -ms + Smim^jC^^Ch (24) 

where Ci = {aj^a^*) is the angular power spectrum. From the above, one also finds that 
the variance of -B3 depends heavily upon the choice of multipole moment pairs {liirii). The 
variance of the terms where h ^ I2 ^ h is given by 

(t| = Ci^Ci^Ci.^, (25) 
9 



while for h = h = I3 = I, rrii = rn2 = ""^s = 0, it is given by 



(T3 — IbCij^Ci^Ci^ — 15Cf. 



(26) 



The variance of the bispectrum can differ by more than one order of magnitude by choosing 
different muhipole pairs. It is easy to see why the cosmic variance of skewness is large: by 
writing it in terms of the angular bispectrum, the variance of the skewness is given by 



This is a sum over all different multipole moments, and is dominated by terms such as 
rrii = OJi = I2 = h = I- The advantage of using the angular bispectrum is that one can 
choose multipole moment pairs to minimize the cosmic variance, and at the same time, to 
maximize model predictions for the angular bispectrum. 

One way to reduce cosmic variance is to subtract out the lower multipole moments. Part 
of the GBR signals will then also be removed so that the signal-noise ratio will be reduced. 
However, this will not be a serious problem for COBE map with fours years of data. The 
noise levels in COBE maps diminish rapidly (oc time^^^^) with additional data (Hinshaw et 
al. 1993). In the four year maps, the noise level will be roughly eight times smaller than in 
the first year skymap. For a scale invariant H-Z initial fluctuation spectrum, the GBR signal 
is reduced by merely a factor of 3 by removing the quadrupole. Thus, the signal-noise ratio 
in the quadrupole-removed four- year map will still be a factor of three better than the first 
year map (with quadrupole). 

A more serious problem with subtracting out lower moments is that part of the non- 
Gaussian signal will also be removed. The amount of non-Gaussian signal removed is model- 
dependent. We list in Table 1 the angular bispectrum B^, modeled by Eq. (|13), in various 
non-Gaussian models, compared with the cosmic variance a^. In non-standard infiation 
or late-time phase transition models, most of the non-Gaussian signal will be removed by 



aliknii, l2m2, hm3)A'^ {likh, 171117121713) . 



(27) 



lirrii 
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subtracting out the lower multipole moments. In these types of model, one should focus on 
the lowest multipole moment (1 — 2). The low hmits on the model constants are: a > 0.6 
in non-standard inflation, where /3 = 7 = 0; /? > 3.2 for late-time phase transitions, where 
a = 0,7 = 1. In models where a — 0,7 = 1/2, the ratio of non-Gaussian signal to the 
cosmic variance is approximately the same for the low / moments, but since the instrumental 
noise is relatively high for large moments, one will also lose information by subtracting out 
lower moments. The low limit on /5 is /3 > 1.5. For most physically motivated models 
(non-standard inflation, non-linear gravitational evolution, the late-time phase transition or 
the 0{N) (7— models) where the model constants a, /3 are less than unity, all terms of the 
angular bispectrum are smaller than the cosmic variance term. Since the cosmic variance 
term oc for the H-Z initial fluctuation spectrum, if the angular bispectrum falls off 
faster than it will be impossible for the experiments to detect deviations from Gaussian 
behavior, unless the amplitude of the non-Gaussian signature is unrealistically large. This 
suggests that the decisive test of the Gaussian nature of fluctuations has to come from degree 
scale experiments, as we have stressed before (Luo 1993). Nevertheless, in models where the 
angular bispectrum is flat (7 = 0), even a very small deviation from Gaussianity (/3 ~ 0.02) 
can stand out above the cosmic variance, by analyzing the angular bispectrum at higher 
multipole moments (Z ~ 8). 

We wish to thank D.N. Schramm, J. Silk, M. White, D. Scott, M. Srednicki and anony- 
mous refree for helpful comments and suggestions. This work was supported in part by DoE 
and NSF. 
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TABLE CAPTION 

Table 1: Comparison of the angular bispectrum ^3 in various non-Gaussian models with the 
cosmic variance of the angular bispectrum for Gaussian fluctuations, weighted by nii). 
The multipole pairs (Zj, mj) are chosen to obtain the maximum value for B3 in non-Gaussian 
models, and at the same time, to keep the cosmic variance as small as possible. 
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Table 1 
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